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CO EMISSION IN OPTICALLY OBSCURED (TYPE-2) QUASARS AT REDSHIFTS z « 0.1-0.4 



<N 

o 

<N 



On 
O 

u 

Oh! 

o 
is 



> 

ON 
i— H 

in 
o 

(N 



X 



M. Krips, R. Neri and P. Cox 

Institut de Radio Astronomie Millimetrique (IRAM), 300, rue de la Piscine, Domaine Universitaire, 38406 Saint Martin d'Heres, France 

Not to appear in Nonlearned J., ^5. 

ABSTRACT 

We present a search for CO emission in a sample of ten type-2 quasar host galaxies with redshifts 
of z^0. 1-0.4. We detect CO(J=1-0) line emission with >5a in the velocity integrated intensity maps 
of five sources. A sixth source shows a tentative detection at the ~4.5<r level of its CO(J=1-0) line 
emission. The CO emission of all six sources is spatially coincident with the position at optical, 
infrared or radio wavelengths. The spectroscopic redshifts derived from the CO(J=1-0) line are very 
close to the photometric ones for all five detections except for the tentative detection for which we 
find a much larger discrepancy. We derive gas masses of ~(2-16) x 1O 9 M0 for the CO emission in the 
six detected sources, while we constrain the gas masses to upper limits of M gas <8x 1O 9 M for the four 
non-detections. These values are of the order or slightly lower than those derived for type-1 quasars. 
The line profiles of the CO(J=1-0) emission are rather narrow (<300 km s _1 ) and single peaked, 
unveiling no typical signatures for current or recent merger activity, and are comparable to that of 
type-1 quasars. However, at least one of the observed sources shows a tidal-tail like emission in the 
optical that is indicative for an on-going or past merging event. 

We also address the problem of detecting spurious ~5<r emission peaks within the field of view. 
Subject headings: galaxies: active — galaxies: ISM — infrared radiation — ISM: molecules — quasars: 
general - radio emission lines 



1. INTRODUCTION 

Molecular gas is an important diagnostic tracer to 
study the activity processes in galaxies, whether in form 
of an accreting super-massive central black hole (AGN) 
or intense starbursts (SB) or both. It does not only build 
up the fuel for the activity but also traces essential physi- 
cal properties of the host galaxy, such as the kinematics, 
distribution, excitation conditions and chemistry of its 
gas. As activity occurs in many different types of galax- 
ies and at various levels of intensity and galaxy evolution, 
a big effort is put into elucidating the physics behind its 
different occurrences. 

Given the elementary importance of molecular gas, 
numerous studies have been conducted to detect and 
map its distribution and kinematics in galaxies. Most 
of them were focussed on sources in the local uni- 
verse given their easier detect ability. However, recent 
technical improvements make it now possi ble to probe 
the hi g her redshift ga l axy populat i on (e.g., IWang et al.l 



(AGN-powered) quasar population still lacks a similarly 
high number of systematic studies. This lack of informa- 
tion is hence a shortcoming in current theories of galaxy 
evolution. 

While most recent molecular gas studies of quasars 
at redshifts of z<0.5 concentrate on opticall y se- 



at reasniits or z<u.o concentrate on opticall y se- 
lected sources (type-1 quasars e.g., IWink et al.l ll997: 
Evans et al Jl200lt iScqville et al.ll200llKrips et al.ll2QQ5a . 



| 2011at iRiechers et al.l I2Q11L 120091: iTacconi et al 
12006 ; Daddi et all 1201(1 120091 : iGreve et al.l 
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iNeri et al.ll2003l ). This is especially important because 
most of the active galaxies, such as the ultra-luminous 
infrared (ULIRGs) /sub-millimeter galaxies and quasars, 
are found in the high-redshift universe emphasizing the 
tight correlation of the activity of a galaxy and its evo- 
lution. However, it is still poorly understood whether 
the nature of the activity plays a role for the evolu- 
tion of a galaxy and how; most current theories as- 
sume that both activity types mark a special stage of 
evolution in a galax y's life (e.g.. [Sanders et aHll988aU bl: 
IHopkins et al.l l2QQ6h . It is thus surprising that molecu- 
lar gas has been extensively studied in a large number of 
(SB-powered) ULIRGs/submillimeter galaxies, while the 
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evidence has recently been 
presented that optically obscured (type-2/reddened type- 
1) quasars may be fou nd at least in eq u al amounts 
(at least for z <4: [R eyes et al. 200c| lLacy et all 
l2QQ7aHMartinez-Sansigre et al.ll2006T ). Identifying type-2 
quasars has been proven to be extremely difficult because 
of their high obscuration at most wavelengths but success 
rates have been significantly improved by several groups 
in the past few years using the less o bscured radio, X- 
ray, and/or Mid-Infra red bands (e.g., ILacv et al.ll2007bl : 
IZakamska er"aTll2008D . 

The study of both, type-1 and type-2 quasars, is criti- 
cal as it challenges the viewing angle and merger-driven 
unification theories for quasars. While in the vie wing- 
angle unification an edge-on clumpy torus blocks the 
direct view onto type-2 quasars (Antonucci 1993), the 
merger-driven unification assumes that type-2 quasars 
are in an ear lier stage of a merg ing event than type-1 
quasars (e.g.. IHopkins et all 12006). The latter scenario 
finds support by recent observational results: 1.) type-2 
quasars are found to have higher median star formation 
luminosities than type-1 quasars suggesting that type-2 
quasar activity goes along with an int ensified star for- 
matio n activity in the host galaxy (e.g., Zaka mska et al.l 
2008) and 2.) the type-2 quasar fraction increases with 
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TABLE 1 
Source list 





Source" 


Redshift b 




Scale 


CU2000 


0J2OOO 


Type" 






2>opt 


(Mpc) 


(kpc/arcsec) 


(hh:mm:ss.s) 


(dd:mm:ss.s) 




SWIRE2 


J021638. 21-042250. 8 


0.304 


1565 


4.5 


02:16:38.20 


-04:22:50.8 


2 


SWIRE2 


J021909.60-052512.9 


0.099 


450 


1.7 


02:19:09.60 


-05:25:12.4 


2 


SWIRE2 


J021939.08-051133.8 


0.150 


706 


2.5 


02:19:39.09 


-05:11:33.8 


2 


SWIRE2 


J022306.74-050529.1 


0.330 


1722 


4.7 


02:23:06.70 


-05:05:29.1 


2 


SWIRE2 


J022508.33-053917.7 


0.293 


1499 


4.3 


02:25:08.30 


-05:39:17.7 


2 


SDSS 


J092014.11+453157.3 


0.403 


2179 


5.4 


09:20:14.10 


+45:31:57.0 


2 


SDSS 


J103951.49+643004.2 


0.402 


2173 


5.4 


10:39:51.50 


+64:30:04.0 


2 


SSTXFLS J171325.1+590531 


0.126 


583 


2.2 


17:13:25.10 


+59:05:31.0 


2 


SSTXFLS J171335.1+584756 


0.133 


619 


2.3 


17:13:35.10 


+58:47:56.0 


1R 


SSTXFLS J172123. 1+601214 


0.325 


1691 


4.8 


17:21:23.10 


+60:12:14.0 


2 



a taken from lLacv et~aTl ([2007b ) and Zakamsk aeFaLI p008h . 
b Optical redshifts (spectroscopically determined). 

c Assuming the following cosmology: Qm=0.27, £1^=0.73 and Ho=71 km s _1 Mpc -1 . 

d Positions taken from lSimpson et al.l pu06h fradioV lLacv et al.l (|2007bl ) (optical/IR) or lZakamska et al.l p008ft 
(optical/IR). 

e Type of AGN: 2= optically obscured AGN; lR=reddened type-1 AGN (type-l= un-obscured AGN). 



in optically obscured (type-2) quasars at redshifts 
0.09<z<0.41. The general properties of the sources in 
the sample are summarized in Table [TJ The sources 
were selected based on their 24/im- fluxes (MIPS, IRAC) 
and spectroscop ic redshifts from th e typ e-2 quasar stud- 
i es con ducted bv lLacv et al.l (|2QQ7bl ) and lZakamska et all 
(2008). The upper limit for the redshifts was chosen so 
that the CO(J=1-0) line emission still falls within the 
3mm band of the IRAM Plateau de Bure Interferom- 
eter (PdBI) rec e ivers. The type-2 quasar sample from 
IZakamska etHI (2008) originates from the SDSS cata- 
logue a nd was followed-up with Spitzer. The two samples 
used in lLacv et al.l (|2007bD are taken from the SWIRE- 
(RA=02h) and the Spitzer XFLS fields (RA=17h). 

2. OBSERVATIONS 

The observations of the CO(J=1-0) emission in a sam- 
ple of ten type-2 quasars (Table [1]) have been conducted 
at the IRAM PdBI in 2007, 2009 and 2010, using five 
to six antennas in the most compact configuration (D) 
with baselines ranging between 15 m and 100 m. A jour- 
nal of the observations is given in Table [2j We tuned the 
3 mm receivers of the PdBI to the redshifted frequency of 
CO(J=1-0) (115.271 GHz at rest) for each source, using 
a bandwidth of 1 GHz (^3000-3500 km/s) and a spectral 
resolution of 2.5 MHz for the observations dated prior to 
2010. In 2010, the new wideband correlator WideX was 
installed at the PdBI so that we had a total of 3.6 GHz 
of bandwidth at our disposal. Typical system tempera- 
tures have been around ~80-200 K during the observa- 
tions. We discarded all data with phase noise exceeding 
50 degrees. We calibrated the data using a set of strong 
sources, including MWC349 for absolute flux calibration. 
The uncertainty of the absolute flux calibration is esti- 
mated to be <10%. 

3. RESULTS 

3.1. Line Emission 

We detect CO(J=1-0) emission in five of the ten type- 
2 quasars and find a tentative detection in a sixth source 
(see Table H and Fig. [Q & EJ, resulting in a -60% de- 
tection rate whi ch is similar to previous CO surveys 
on quasars (e.g., lEvans et al.l 120011 : iScoville et all 120031 : 



iBertram et alll2007h . The CO(J=1-0) emission in the 
remaining four sources remains undetected. 

In the course of our observations, we "detected" >5a 
emission peaks at off-center positions in the field of view 
of at least four of our ten targets. None of these peaks 
were found to coincide in position with any of our quasars 
nor any known galaxy in the NED data base. Statistical 
considerations (see Appendix) further argue for a spu- 
rious nature of these ~ 5<r peaks; we dedicate a section 
on the Appendix to this problem. Otherwise no further 
consideration is given to these noise peaks in the follow- 
ing. 

3.1.1. Detections at the Phase Center 

The line emission of five of our ten targets is clearly 
detected in the line spectrum and the integrated maps. 
The spatial position of the emission (Fig. [2]) as well as 
the redshift derived from the line centers (Fig. [I]) coincide 
well with t he spatial positions fr om optical/IR/radio ob- 
servations ([Simpson et al.ll2006l Table [3]) and their pho- 
tometric redshifts (see Table [1] and [4]) . The CO redshifts 
were determined by fitting a Gaussian profile to the line 
spectrum of each source. 

We further find a tentative detection of the CO(J=l- 
0) emission in a sixth source (see Fig. [T]and[2j). A 4.5<r 
peak is found close to the position of the radio emission 
of this source. The CO line, however, is significantly 
shifted bluewards with respect to what is expected from 
the photometric redshift. The difference is much larger 
than for any of the other five detected sources. Also, 
the line profile of the CO(J=1-0) emission is not statis- 
tically significant, so that additional observations would 
be needed to confirm or discard this tentative detection. 
However, since the 4.5<r peak is within the radio position 
of this quasar, we think that it is reasonable to assume 
that this detection is indeed real but we will still treat 
this source as tentative for the remainder of the discus- 
sions. 

All six sources exhibit simple line profiles that appear 
Gaussian in nature, single peaked and narrow. To derive 
the line intensities, line centers, and line width, we fitted 
a Gaussian function to all six lines. The results are given 
in Table gj 

Four out of these six galaxies appear to be unresolved 
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Fig. 1. — Spatially integrated CO(J=1-0) line spectrum of all six CO detections. The zero velocity corresponds to the (spectroscopically 
derived) optical redshift (dotted grey lines). The new redshifts based on the CO(J=1-0) lines are determined from fitting a Gaussian profile 
to the spectra (dashed black lines). 



TABLE 2 
Journal of Observations 



Source 


Observing Date 
(Year: Month) 


Obs. Frequency 
(GHz) 


Integration time a 
(hours) 


Synthesized Beam b 
(" x"@°) 


J021638-042250 


2009:May,Jun,Jul 


88.398 


7.4 


7.4 x 5.0(Q 


6 153 


J021909-052512 


2007:Oct,Nov 


104.887 


1.9 


6.9 x 4.8(Q 


D 21 


J021939-051133 


2007:Sep 


100.236 


4.0 


6.0 x 4.7(Q 


D161 


J022306-050529 


2009:May,Jun 


86.670 


10.3 


8.3 x 5.1 <§ 


6 148 


J022508-053917 


2009:Sep,Oct 


89.150 


5.4 


11.1 x 4.6(6 


^146 


J092014+453157 


2009:May,Jun,Aug; 2010:May,Apr 


82.206 


11.3 


6.2 x 4.7 <§ 


6 127 


J103951+643004 


2009:May,Jun; 2010:Apr 


82.219 


7.6 


6.5 x 5.4 @ 


6129 


J171325+590531 


2007:Aug 


102.372 


3.0 


5.9 x 3.8 (6 


6 74 


J171335+584756 


2007:Oct; 2009:Sep 


101.740 


8.2 


5.2 x 3.3 <§ 


6 69 


J172123+601214 


2007:Aug,Sep,Oct 


86.997 


6.5 


7.2 x 5.2 (6 


6 134 



a On-source integration time estimated with respect to 6 antennas and dual polarisation. 
b For natural weighting. 



at their detection level and no significant velocity gra- 
dient can be identified. Only the CO emission regions 
of J021909-052513 and J021939-051133 appear to be 
slightly extended and show a velocity gradient along a 
North-East to South- West direction (Fig. [3]) although the 
one for J021909— 052513 is less pronounced. W hile the 
optical i'-image (taken from lSimpson et aD2006[ ) shows a 
tail-like extension for J021909— 052513 indicating a possi- 
ble interaction with another galaxy, the optical image for 
J021939— 051133 shows a compact galaxy with no signs 
of disruption. The CO emission in J021909— 052513 sug- 
gests the presence of molecular gas along the tidal tail 
of this merger system. Fitting a Gaussian to the uv- 
data results in a size estimate of 2±1" (—5 kpc), i.e., 
very centralised molecular gas. However, at the angular 
resolution and sensitivity of our observations, we cannot 
exclude that part of the more extended diffuse gas might 
have been resolved out. 

3.1.2. Non- detections 

Although we cannot entirely rule out this possibil- 
ity, we doubt that the four non-detections are due 



to uncertainties in t he optical spectroscopic redshifts. 
lHainline et al.l (|2004t ) present a comparison between op- 
tical and CO redshifts in high-redshift sources . The 
authors find, excluding a few outliers, a difference of 
^ z max = | z opt~ zco 1^0-008, which corresponds to a ve- 
locity shift of -±(1500-2000) km s _1 at 3.5 mm in the 
redshift range z=0.2-0.5. A similar deviation of —0.008 
is seen wh en adding the CO detections at high -z since 
2005 (e.g., iGreve et aLll2QQ5t ICoppin et al.ll2008[ ). Inter- 
estingly, the outliers which show 0.008<Az<0.05 all lie 
at redshift z>2.5. Therefore, the large bandwidth of the 
observations of 3000 km s _1 should be sufficiently large 
to take out uncertainties in the optical redshift for the 
sources presented in this paper. However, it has to be 
pointed out that for the determination of the difference 
Az only sources with detected CO emission were consid- 
ered. The actual deviation Az could be hence larger than 
assumed here. Only larger bandwidth (>4 GHz) and 
higher sensitivity (going down to M gas <lO 8 M0 where 
M gas =M(H2+He) is the total gas mass) observations will 
allow us to reliably eliminate all uncertainties from opti- 
cal redshifts for the detection of CO emission in higher-z 
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Fig. 2. — Velocity integrated CO(J=1-0) line emission maps of the six detected type-2 quasar host galaxies. The large dashed circles 
indicate the field of view for each target. The synthesized beams are shown in the lower right corner of each map. Negative contours 
(dotted) correspond to — 3,— 2, — lcr, and positive contours (solid) start at la in steps of la with the following noise levels for each map: (a) 
lcr^O.05 Jy beam" 1 km s _1 (Av int ~220 km s _1 ); (b) lcr=0.58 Jy beam" 1 km s _1 (Av int ~490 km s _1 ); (c) lcr=0.39 Jy beam" 1 km s _1 
(Av int ~330 km s" 1 ); (d) lcr=0.08 Jy beam" 1 km s" 1 (Av int ~400 km s" 1 ); (e) lcr=0.08 Jy beam" 1 km s" 1 (Av int ~350 km s" 1 ); (f) 
lcr^O.07 Jy beam" 1 km s _1 (Av int ~420 km s _1 ) 



3.2. Continuum Emission 

Continuum emission at 3 mm has not been found in 
any of the ten sources. The 3<r upper limits derived 
for the 3 mm continuum fluxes (Table |4]) are similar or 
slightly lower than the 1.4 GHz flux densities (Table [3]) 
excluding a flat radio spectrum in these sources and 
pointing rather towards a steep radio spectrum. A steep 
radio spectrum is generally associated with non-thermal 
emission which is typically found in this galaxy popula- 
tion. However, extrapolating the 1.4 GHz flux densities 
(Table [3]) with a spectral index of —1 (i.e., f^oc^ -1 ) in- 
fers flux densities at 3 mm that are smaller by an order 
of magnitude than our 3 sigma upper limits. We also 
stacked the continuum maps from all ten targets with no 
significant detection in the stacked image, i.e., the emis- 
sion peak in the center remains at <3<t=0.09 mJy. Thus, 
no definite conclusion can be drawn on the nature of the 
3 mm continuum emission in our sources. 

4. DISCUSSION AND CONCLUSIONS 

4.1. Molecular Gas Masses 

Using the standard Galactic conversion factor for the 
Mgas-to-L'co ratio of X C o=4.8 M (K km s _1 pc 2 ) -1 



fe.g.. lSolomon fc Barretl fl991). we find a range of the to- 
tal gas mass of M gas ~(2-16) x 1O 9 M for all five detected 
sources (plus the one tentatively detected) and 3<r upper 
limits for the four non- detections of M gas <(2-8) x 1O 9 M 
(Table 2]). However, the standard Galactic conversion 
factor Xco is known to over-estimate the total gas mass 
in the ULIRG / su bmillimeter population b y more than a 
factor of 5 (e.g.. iDownes fc Solomon! fl998[ ) and seems to 
be also highl y dependent on t he physical properties (e.g., 
metallicity; IBell et al.ll2007all ). iBell et al.l (|2007bh find 
indications for a lower Xco in AGN dominated galaxies 
in the local universe as well (e.g., M51). If this were also 
true for higher-z quasars, the gas masses for the sources 
in this paper would reduce to values below 3xlO 9 M . 

4.2. Dynamical Masses 

Assuming a radius of the emission equal or lower than 
half the synthesized beam (i.e., ~3-5"~5-15 kpc) and 
taking the full width at half maximum (FWHM) of the 
line profiles determined for the six detected sources, the 
dynamical masses amount to ~(1O 1O -1O 12 )M xsin _2 z 
with i being the disk inclination angle. These dynami- 
cal masses easily account for the obtained total gas mass 
within the assumed radius of 5-15 kpc (~ 10% of the 
stellar amount in the entire galaxy, which is a few times 
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TABLE 3 

BOLOMETRIC (RADIO, IR, X-RAY) PROPERTIES OF THE OBSERVED SOURCES 



Source 


Sl60/Lim a 


S70/j,m a 


S24/^m a 




Lfir b 


Lmir b 


Lir b 


Sl.4GHz C 


So. 2-12 keV d 




(mJy) 


(mJy) 


(mJy) 


(mJy) 


(1O 1O L ) 


(1O 1O L ) 


(1O 1O L ) 


(mJy) 


(10- 14 ergs" 1 cm- 2 ) 


J021638-042251 




27 


14.6 


2.7 


>8 


-23 


>32 






J021909-052513 


152 


130 


25.6 


3.1 


-6 


-3 


-9 


1.36±0.02 


2.6 


J021939-051133 




70 


32.9 


8.5 


~4 


-11 


>15 


1.24±0.06 


5.3 


J022306-050529 






15.7 


1.5 




-30 


>30 




0.95 


J022508-053917 




32 


9.6 


2.6 


>9 


-14 


>23 






J092014+453157 














-8 




44 


J103951+643004 














-13 


4.3±0.5 


70 


J171325+590531 






9.5 


1.19 




-2 


>2 


0.16±0.03 




J171335+584756 


146 


97 


23.7 




-9 


-6 


-15 


3.4±0.2 




J172123+601214 






13.3 


3.71 


~4 d 


-25 


-29 


0.26±0.04 





a taken from ILacv et al. (2004), Lacv et aTl (l2007bT L [Fraver et al. (120061 ) and/or from the SWIRE catalogue available at 
http : / / swire . ipac . caltech . edu/ swire/ astronomers/ data_access . html . 



b dete rmined based on LFiR=C2(z)^L^(70/im)+C3(z)^W(160^m), L]y[iR=Ci( 2; ) z/ L Z y (25/im), and Lir=Lfir+Lmir from Dale Sz Heloul 
(US!); with &(z « 0.1)^1.75, Ci(z w 0.3)^2.2, ^{z w 0.1) » Q 2 (z ^ 0.3)^0 82, Cs jz - 0.1)«C 3 (s - 0.3)^1.347. The IR luminosities 
for J092014+453157 and J 103951+64300 4 were taken from Zakamska et al. (120081). 

c References: J17-sources: ICondon et"aH (|2003T ), J02-sources: ISimpson et al.l ((2006), J103951+643004 from lCondon et al.l (1998). 
d from ILacv eTaTl (l2007bT ) and Zakamsk a"eTaII (120081 ) . 



TABLE 4 
Observational Results 



Source 


RMS a 


f b 
icont 


fco c ' d 




zco 9 


Sco d 


L'co d 


M gaB d > h 




(mJy) 


(mJy) 


(mJy) 


(km s - 1 ) (km s - 1 ) 




(Jy km s 1 ) 


(K km s 1 pc 2 ) 


(1O 9 M ) 










DETECTIONS 










J021638-042251 


0.8 


<0.3 


2.0±0.4 


170±40 +102±17 


0.3044±0.0002 


0.5±0.1 


(23±5)xl0 8 


11±2 


J021909-052513 


1.4 


<0.6 


5.6±0.6 


300±40 -135±16 


0.0985±0.0006 


1.5±0.3 


( 7±l)xl0 8 


3.2±0.6 


J021939-051133 


1.6 


<0.5 


7.3±0.7 


200±30 +100±10 


0.1504±0.0002 


1.6±0.2 


(17±2)xl0 8 


8±1 


J022306-050529 


0.7 


<0.2 


1.7±0.3 


290±60 +98±27 


0.3304±0.0004 


0.6±0.1 


(33±5)xl0 8 


16±3 


J171335+584756 


0.8 


<0.3 


2.5±0.4 


190±30 +480±10 


0.1348±0.0003 


0.45±0.08 


( 5±l)xl0 8 


2.4±0.4 










TENTATIVE DETECTIONS 








J103951+643004 


0.6 


<0.2 


0.81±0.26 


310±120 -1415±50 


0.3954±0.0006 


0.25±0.06 


(21±5)xl0 8 


10±2 










NON-DETECTIONS 








J022508-053917 


0.3 


<0.3 


<0.9 






<0.3* 


<13xl0 8 


<6 


J092014+453157 


0.2 


<0.2 


<0.6 






<0.2* 


<17xl0 8 


<8 


J171325+590531 


0.5 


<0.5 


<1.5 






<0.4* 


<3xl0 8 


<2 


J172123+601214 


0.2 


<0.2 


<0.6 






<0.2* 


<llxl0 8 


<5 





















a RMS noise determined from channels with 10 MHz spectral resolution (~28 km s 1 ) for the CO detections and with 105 MHz spectral resolution 
(~300 km s" 1 ) for the CO non-dectections. 

k 3cr-upper limits for the continuum flux averaged over the line-free channels for each source. 
c Peak line fluxes for the detections and 3<r-upper limits for the non-detections. 

d Errors are of pure statistical nature and do not account for uncertainties from the absolute flux calibration. The latter is estimated to be at a 
conservative level of ~10%. The integrated flux densities have been determined assuming that the continuum emission in these sources is negligible 
at the given noise level (see also text). 

e Avfwhm= Full Width at Half Maximum of the line. Determined from CO line profile by fitting a Gaussian to the data. 
^ Systemic velocity; zero velocity corresponding to the optical redshifts of the sources. 
9 Redshift determined from the central velocity of the CO(J=l — 0) line. 

The standard galactic M gas -to-L'co conversion factor of X C o=4.8M (K km s _1 pc 2 ) -1 has been assumed. 
1 Assuming a Ai>fwhm of ~260 km s _ 1 corresponding to the average of the FWHM of the detected sources. 



1O 12 M0) and leave enough room for additional mate- 
rial, such as the super-massive central black hole (up to 
~1O 1O M0 in the most extreme cases) and the stellar con- 
tent (up to a few times IO^Mq). 

4.3. CO versus IR Luminosities 

Comparing the CO-luminosities for the detected 
sources and the upper limits for the undetected ones with 
those of type-1 qu a sars taken from the l iterature (e.g., 
lEvans et al.l 1200 It iScoville et all 120031 : iBertram et all 
120071 ). we find that they fall within the same area in 



the L'co-Lir diagram (Fig. 3]). However, it has to be 
noted that most of the IR luminosites assumed for the 
type-2(lR) quasars must probably be regarded as lower 
limits (see Table [3]). Several of our type-2 quasars have 
no measured/detected FIR fluxes. Also, obscuration 
through dust might significantly reduce the emission at 
MIR wavelengths. However, based on comparison with 
other dust enshrouded quasars, the IR-luminosities can 
be larger by up to an order of magnitude than the MIR- 
luminosities reducing hence the significance of the MIR 
to the total IR- luminosities. 
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Fig. 3. Blue- and redshifted CO(J=1-0) line emission (contours) of J021909-052513 (left panel) and J0219 39-051133 (rigt panel ) 
overlaid to the Subaru/Suprime-Cam i' image (grey scale) and the radio emission (faint grey contours) taken from Simpson et al. (2006 
The velocity gradient is identicated with a green line; the green line is dashed for J021909— 052513 because a velocity gradient is much 
more difficult to identify than for J021939— 051133. Contours start at 2a and go in steps of 2cr. 




Taking all published CO detections in quasars into 
account and excluding the upper and lower limits for 
some of the CO and IR luminosities, we find a lin- 
ear correlation between the CO and IR luminosities 
of 0.70±0.03 dex (Fig. 2J. When considering the up- 
per limits for the CO luminosities and, separately, the 
lower limits for the IR luminosities, we obtain a very 
similar slope of 0.69 ±0.03 dex based on t he an alytical 
method described in Fei gelson fc Nelsonl (|l985f ). This 
is in good agreement with resul ts found for ULIRGs 
and nearby active galax ies (e.g.. iGao fc Solomon! l2Q04t 
iNaravanan et al.l l2008al ). However, some of the upper 
limits of the CO observations fall significantly below 
the correlation but it is unclear whether this is due 
to a breakdown in the correlation, to uncertainties in 
the optical redshift and/or the significant contribution 
of the AGN to the (M)IR-luminosities that might even 
dominate over that from s tar-formation (e.g. jLacv et al.l 
2QQ7bl: [Brand et all [2006). Besides obscuration effects, 
the latter represents an important factor in determin- 
ing "accurate" (i.e., purely star- formation related) IR- 
luminosities in particular for type-2 quasars. However, 
the significance of the AGN contribution to the energy 
output remains also strongly debated for UL IRGs (e.g., 
iFranceschini et al . 2003; I manishi et al.ll2007l ) so that the 
importance of the AGN influence on the correlation be- 
tween CO and IR-luminosities is hard to quantify. On 
the other hand, the similarity of this correlation between 
ULIRGs and quasars may indicate that either the AGN 
contribution influences the correlation in the same way 
for these sources (such as adding an offset in a system- 
atic way) or is even a fundamental part of the correlation. 
A detailed entangling of the different components is cer- 
tainly necessary to further our understanding of this very 
fundamental correlation. 

Alternatively, the slope of the L'co-Lir correla- 



tion coul d also vary for different galaxies (as dis - 
cussed bv IGao fc Solomon! I2004J : iBertram et all l2QQ7h . 
iBertram et al.l (|2007f ) suggests a scenario in which ge- 
ometrical effects naturally lead to a change in slope: 
galaxies with highly centralised and confined molecular 
gas (such as ULIRGs) have a larger surface filling factor 
and thus a steeper slope whilst 'normal' weakly interact- 
ing/active galaxies show more extended molecular gas 
resulting in a milder slope. The similarity of the slope 
between ULIRGs and quasars might hence be an indica- 
tion that merger events are also important for the quasar 
population. However, more observations and studies of 
the distribution of CO in quasars, and active galaxies in 
general, have to be undertaken to validate this interpre- 
tation. 

4.4. CO Luminosities versus redshift 

Plotting the CO luminosities of the quasars detected 
in CO against their redshifts (right panel of Fig. [4j), sug- 
gests that, despite similar IR luminosities, type-2 quasars 
have lower CO luminosities, and hence a lower amount 
of molecular gas at comparable stages of evolution than 
type-1 quasars, i.e., similar redshifts; this impression is 
mainly given when looking at redshifts between 0.08 and 
0.2 (above redshifts of 0.2 and below redshifts of 1.0 
no CO detection for type-1 quasars has yet been re- 
ported). When applying the stastistical analysis from 
iFeigelson fc Nelsonl (jl985f ) that includes the upper limits 
(see also previous Section) we do find marginally different 
mean values for the CO luminosities in type-1 quasars of 
<L'co>=(i3±4)xl0 8 (K km/s pc 2 ) and those in type-2 
quasars of <L' C o> = (8±3) x 10 8 (K km/s pc 2 ) for red- 
shifts between 0.08 and 0.2; for the entire type-1 and 
type-2 quasar samples we find almost the same mean 
values of <L' CO >=(12±3) x 10 8 (K km/s pc 2 ). How- 
ever, not only do th e mean values agree within the errors, 
a logrank analysis ( Feigel son fc Nelsonl 1 1 98 b\ ) also shows 
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Fig. 4. — CO versus IR luminosity plot (left) and CO luminosity versus redshift plot (right) of quasars with CO measurements taken from 
this paper (filled black circles) and the literature (open symbols). The solid black curve (left) is a linear fit to the data. Given their high 
dust obscuration, the IR-luminosities of the type-2(lR) quasars can only be regarded as lower limits. Upper and lower limits are marked 
with arrows. For some of the hig h-redshift source s (z>l), only FIR luminosities were available, so we assume Ljr=1.5xLfir based on a 
LnR-to-LiR comparison done bv lPott e t al. (2006J); please note that some quasars can show a ratio larger by a factor of 2-4 because of the 
contribution of the AGN to the MIR luminosity. However, increasing the ratio to 4 does not significantly change the slo pe of the cor r elatio n 
and is taken into ac count in the err or of the fit. The refe rences for the litera ture data ar e as foll ow s: PG quasars — Evans et al. (2001); 
Scoville et al. (2003); HE quasars - Bertram et al. (2007); high-z quasars - Wang et al. (2011a. b); Polletta et al. (2011); Coppin et al. 
(2008); Aravena et al. (2008); Maiolino et al. (2007); Krips et al. (2005b); Walter et al. (2004); Carilli et al. (2002); Cox et al. (2002). 



that both samples are similar at a 50% significance level. 
Therefore we can neither rule out that both samples ex- 
hibit similar CO luminosities nor can we support it with 
the current samples. 

Even if the CO luminosities were proven to be lower 
in type-2 than type-1 sources this would be at odds with 
both unification theories. The nature of the viewing- 
angle unification theory suggests that no statistically rel- 
evant differences should exist between CO luminosities 
of type-1 and type-2 quasars. On the other hand, if the 
merger-driven unification were to hold true, higher CO 
luminosities (i.e., higher molecular gas masses) would be 
expected for type-2 quasars. 

Thus, different explanations must be considered. One 
would be the low number statistics. An alternative 
reason could be the different selection criteria used for 
the CO surveys of the type-1 and type-2 quasars. It 
seems that the type-1 quasars were chosen based on their 
bright IR luminosities while this wasn't a criterium for 
the type-2 sources. The type-1 quasar sample might 
hence be biased toward gas rich sources. Also, we can- 
not exclude resolution effects and hence an underesti- 
mation of the CO luminosities in our interferometric 
study. Ho wever, the interfer omet ric observations car - 
ried out bv lEvans etHI (|2QQlh and iScoville et al.l (|2QQ3h 
probably suffered from similar problems given the similar 
redshift ranges of the targets and hence similar spatial 
scales. Only larger samples will help to settle the ques- 
tion whether the CO luminosities are different for the 
two quasar populations. 



4.5. Linewidth of CO Emission 



In case of a lack of sufficient angular resolution, 
the linewidth of molecular lines can be a first-order 
merger-activity indicator for galaxies, i.e., the larger the 
linewidth the more likely a galaxy might be interact- 
ing or merging with a nother galaxy (e.g., Arp 220; see 
also lGreve et al.l l2005): however, the contrary is not nec- 
essarily true. Comparing hence the linewidths of our 
type-2 quasars with that of ULIRGs and high-z galaxies 
(Fig. [5]), for which a high percentage is assumed to be in 
a merger stage, shows indeed that the distribution of the 
linewidth peaks at higher values for the ULIRG/high- 
redshift galaxy population than for the nearby type-2 
quasars. Also, the ULIRG/high-redshift galaxy popu- 
lation exhibits some very large linewidths with values 
exceeding 500 km/s as opposed to the type-2 quasars for 
which all measured linewidths stay well below 400 km/s. 
No clear difference though can be identified between 
the type-2 and type-1 quasars; the same is also true 
for the more local type-1 and type-2 Seyfert galaxies. 
Both quasars and Seyferts reveal a similar distribution 
of their linewidths. A significant difference in the distri- 
bution of the linewidth between the type-1 and type-2 
quasars could favor a merger-driven unification theory; 
the broadest linewidths should occur for pre-coalescent 
mergers for which the gas disks of the merging galaxies 
have not yet met (e.g. Jlono et al.lfeOOQh . 

Even in the viewing- angle unification, one might ex- 
pect a difference between the linewidths for type-1 and 
type-2 quasars under certain assumptions. Let us as- 
sume that the molecular gas is mostly located in a disk 
in these quasars and that this (large-scale) disk exhibits 
a similar inclination to the putative (small-scale) torus. 
If that were true, the CO emission would be mostly seen 
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Fig. 5. — A comparison of the measured CO linewidth of the 
six detected type-2 quasars to other active galaxies in the local 
and high-redshift universe; the histogramms were derived such as 
that FWHM min < FWHM < FWHM max . The FWHM for the 
ULIRGs and high-redsh i ft gal axies are taken from references given 
in Fig.|4]and Iono et al. (2009) , for th e nearb y Seyfert galaxies an d 
type-1 quasars from Barvainis et al. (1989); Alloin et al. (1992); 
Vila-Vilaro et al. (1998); Sakamoto et al. (1999); Evans et al. 
(2001); Scoville et al. (2003); Krips et al. (2005a); Bertram et al. 
(2007). Note of caution: not corrected for inclination; the CO for 
the high-z sources is mostly observed at J U pper >1 so that their 
FWHM might be even higher. 

edge-on in type-2 and face-on in type-1 quasars. In turn 
this means that the CO linewidth would generally be 
larger for type-2 than for type-1 quasars. However, over- 
all the linewidths are not expected to exceed values of 
500 km s _1 as found in mergers. Therefore, a distinction 



to the merger unification should still be possible. On 
the other hand, the assumptions made above may not 
be valid. If a merging event is indeed taking place in 
these systems, the molecular gas might be located in a 
severely disrupted disk for which an inclination would be 
difficult to derive. Also, even if located in a well defined 
(large scale) disk, its inclination does not have to be the 
same as that of a small scale torus (se e for instance the 
discussion for M51 and NGC 1068 in: iMatsushita et al.l 
120071 : iKrips et al.ll201l[ ). 

However, one has to keep in mind that there are a few 
caveats to be considered for this comparison. First, the 
size of the sample suffers from small number statistics. 
Second, no correction has been done for disk inclination 
effects, but, since this is true for all galaxy populations 
used here, the effects are assumed to be on average the 
same for all and to cancel each other out. Third, the CO 
lines for the high-redshift galaxy population have been 
detected in higher rotational transitions which can un- 
derestimate the actual linewidth as they tend to trace the 
warmer and more compact molecular gas. Fourth, also 
the possible presence of gas outflo ws is to be considered a 
linewid t h broadening; factor f e.g.. [Narayanan et al.ll2 006. 

iGreve fc Sommer-Larsenl 120081 : iFeruglio et al.l 



5. SUMMARY 

We detected the CO(J=1-0) line emission in five out 
of a sample of ten type-2 quasars. We further report on 
the tentative detection of the CO(J=1-0) line emission in 
a sixth source. The estimated molecular gas masses are 
rather moderate with values not exceeding 2xlO lo M . 
All detections nicely follow the L'co-to-LiR correlation 
for quasars. The CO luminosities derived for the six de- 
tections and the upper limits for the remaining four type- 
2 quasars seem to lie slightly below that of previously 
detected type-1 quasars for similar redshifts. However, 
given the low number statistics of the current samples 
this difference is not found to be at a statistically signif- 
icant level and might simply be due to a selection bias. 
The rather simple CO(J=1-0) line profiles and moderate 
CO luminosities do not seem to support a recent or cur- 
rent merger activity in these sources; however, at least 
in the case of one source, J021909— 052513, the presence 
of a tidal tail in the optical images is strongly suggestive 
of a recent or on-going merger event. The data discussed 
in this paper do not allow to favor either the merger- 
driven or the viewing angle unification models. Higher 
number statistics are mandatory to support one or the 
other scenario. 



APPENDIX 

DETECTIONS OF ~5<j PEAKS AT OFF-CENTER POSITIONS 

The CO detections and non-detections for the ten sources were determined by running an automated procedure on 
the calibrated data cubes in both the uv- and the (dirty) image plane. This procedure generates integrated maps for 
all consecutive channel combinations, i.e., for a 100-channel data cube the code integrates channel number 1 to 2, 1 
to 3, 1 to 100, 2 to 3, 2 to 4, 2 to 100, until it reaches 99 to 100. It then searches for the emission maximum in 
each of the integrated maps and derives the absolute maximum of those individual maxima. This provides a robust 
and unbiased approach in the search for line emission in large spectral data cubes. With this method several >5<r line 
emission peaks were observed in addition to the detections and non-detections at the center of the arrays field of view. 

The uncertainties of the optical/IR positions of the type-2 quasars are smaller than 1" . This is much smaller than 
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Fig. 6. — Statistical analysis of the Gaussian noise peak distribution for a uniformly filled uv-coverage with one visibility of constant 
weight every 15m. 100 independent simulations for a 100-channel data cube have been run (for an observing frequency of ~100 GHz in D 
configuration of the PdBI). The maps show the average number count per pixel for each peak intensity normalised by the beamsize within 
a certain cr-range. The sum of the total number of peaks per map divided by the number of repetitions and normalised by the beamsize 
is given in the upper left corner of each map. The top panels show the positive peaks, the bottom panels the negative ones. The circles 
indicate radii starting at 0.5 x0 s , S and then increasing in steps of S (0 s =S.5 ,f ). The synthesised beam 6 S is shown in the lower right 
corner. 

the difference between the actually targeted type-2 quasars and the detected emission peaks which were found to be off 
by more than from the quasar position (>2-lO0 s ). A positional search close to the peaks using NED did not show 
any other galaxies that could be associated with these peaks raising doubts about the reliability of these detections. 

Although a >5a peak in the velocity integrated intensity map is generally assumed to mark a reliable detection limit 
for line emission with the IRAM PdBI, even for off center positions, the significant increase in spectral bandwidth 
might necessitate a refinement of this detection limit based on statistical considerations. The probability of finding 
several consecutive 2-3a noise peaks in adjacent spectral channels at any position within (1.5 x) the primary beam 
increases with the spectral bandwidth and ultimately leads to spurious detections. In order to assess the probability 
of spurious detections, we conducted a statistical analysis of the data based on an ideal noise distribution as a starting 
point. Ideal in this context stands for a uniformly filled uv-plane (i.e., one visibility with a constant weight every 
15m in either direction up to the maximum D-configuration baseline of 100m used in the actual observations) and a 
Gaussian noise distribution. We used 100 channels per simulated data cube with 128x128 pixels per channel map. 
We ran the aforementioned "peak-finder" procedure and repeated this 100 times, every time with a re-computed noise 
distribution. The results are shown in Fig. [6] A peak-count map is shown for different levels of detection in units of 
a, i.e., as multiples of the rms noise (starting from 0.5a up to 6.5a peaks). The upper panel shows the average counts 
of positive peaks per pixel normalised to the beamsize while the lower panel that of the negative ones. The absolute 
average count that is derived from the total counts per map divided by the number of iterations (100) and normalised 
by the beamsize (~30 pixel 2 per beam) is plotted in each SNR range map. The circles indicate different radii, starting 
at r=#2/2,#2 and increasing then in steps of S ; in our example the synthesized beam has a size of S =3.5" (for an 
observing frequency of ~100 GHz in D configuration). 

Even in the ideal case of uncorrelated noise, constant visibility weights and even sampling of the uv-plane, we obtain, 
on average, at least one ~5a peak anywhere in (1.5 x) the primary beam. Even the probability to find a peak with 
>5.5a is non-zero and lies at the ~9% level. In other words, for a CO line survey of 10 sources this hence means that 
one spurious >5.5a peak will appear within the field of view of at least one of the ten sources. Furthermore, close 
to the center of the field of view, i.e., within a radius equal to the size of S , the probability to find a ~5a peak is 
estimated to be at the ~5% significance level. Decreasing this circle however further below a radius of ~1.75" (=5 
times the positional uncertainty a p for a synthesized beam of ^ s =3.5 // and a signal to noise ratio of SNR=5a assuming 
c"p=^s/(2SNR)), the probablility drops to almost zero to detect a noise peak with >4.5a, giving statistical support for 
the tentative detection of the CO(J=1-0) emission in J103951+643004. 

Obviously, an uneven sampling of the uv-plane, systematic calibration uncertainties, and correlated noise between 
adjacent channel maps and pixels is likely to increase the number of >4.5a peaks in the field of view. Spurious 
detections will certainly challenge future detection surveys as bandwidths will steadily increase in the coming years, 
especially in the context of a blind survey for line emission from galaxies at high redshift. Based on these statistics, 
we repeated the observations for those sources for which we found >5a off-center peaks within the field of view. None 
of these peaks could be confirmed to be real. 
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